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^ I The time evolution of correlation functions in statistical systems is de- 

scribed by an exact functional differential equation for the corresponding gen- 
erating functionals. This allows for a systematic discussion of non-equilibrium 
physics and the approach to equilibrium without the need of solving the non- 
linear microscopic equations of motion or computing the time dependence of 
q ■ the probability distribution explicitly. 
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A description of the time evolution of statistical systems is needed in very diverse 
areas of science, ranging from the growth of density fluctuations in the early universe 
to the evolution of a population of bacteria. As a typical problem we may consider 
a given distribution of density fluctuations in the early universe at the time when 
the photons decouple - say a Gaussian distribution around central values given by a 
flat Harrison- Zeldovich spectrum - and ask how it looks like when the galaxies form. 
For small fluctuations the underlying gravitational equations of motion (or field 
equations) can be linearized and approximate solutions found |TJ. For large density 
contrasts, however, the nonlinearities in the field equations become important and 
the analytical discussion becomes more involved. Most approaches to this problem 
therefore use numerical simulations. The observed galaxy correlation function shows 
a power-like decay and similar for clusters. This raises the challenge to understand 
the exponents analytically - perhaps by the attraction of the nonlinear system to a 
suitable partial fixed point, similar to the renormalization flow of couplings in the 
theory of critical phenomena. 

Examples from particle physics include the rate of baryon number- violating pro- 
cesses at high temperature and the discussion of a possible formation of a disordered 
chiral condensate in heavy ion collisions. After such a collision the distribution of 
pion fields is not given by the equilibrium distribution corresponding to the ground 
state. The time evolution of the non-equilibrium distribution could produce spec- 
tacular coherent effects. In statistical physics one may ask how a system evolves 
to equilibrium after a quench or starting in a metastable state. We propose here 
a method how such questions can be addressed directly on the level of correlation 
functions without solving the (nonlinear) equations of motions or computing the 
time evolution of the probability distribution. It is based on an exact functional 
differential equation for the time evolution of generating functionals for the correla- 
tions functions. The microscopic laws or equations of motion reflect themselves in 
the precise form of this evolution equation. Many general aspects, as the approach 
to possible fixed points for t — > oo, are, however, independent of this precise form. 

Consider a system of degrees of freedom Xmif) whose time evolution is determined 
by a differential equation of motion 

~nXm = Xm — Fm[x] = 5^ fmn 1 ...n k Xn 1 ---Xn k (1) 

ai k=0 

(We always sum over repeated indices and use X-n — X* n with a summation over 
positive and negative values of n if Xn is a complex variable. In case of continuous 
field equations for xi x i t) the index summation corresponds to an integration over 
x or a corresponding momentum summation for Fourier modes. The equation of 
motion is nonlinear unless vanishes for k > 1 and F m can also depend on 
t.) We assume that with given initial values x m = Xm(to) ec L- (1) has a solution 
Xm{t', X°)i but we will not need its explicit form. The initial probability distribution 
~ exp — So[x°] at time to describes an ensemble of initial conditions. The correlation 
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functions for the variables Xm at some time t > to are then given by 

< Xn^-XnM >= Z' 1 J D X ° XnA^X^-Xn^X^ew-Soix } (2) 

with Z = J Dx° exp —Sq[x°} and / Dx° = fUm^Xm- Inserting the equation of mo- 
tion we see that the evolution of the mean value of Xm involves the higher correlation 
functions 
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< Xm(t) >= E f£L.n k < XnAt)-Xn k (t) > (3) 
Ql k=0 

Writing down similar equations for the correlation functions one arrives at an infinite 
system of differential equations for the time dependence of the n-point functions. 

Insight in the behaviour of this system can be gained from the generating func- 
tional for time-dependent n-point functions 

Z\j,t] = J Dx°exp{-So[x ] +f m X m (t;X°)} (4) 

or the related functionals for the connected correlation functions W\j, t] = In Z[j, t] 
or the 1PI irreducible correlation functions 

V[(p,t] = - In Z[j,t] + j* m p m , 

V™ = o-» (5) 

We will concentrate here on the time evolution of the non-equilibrium effective 
action T[tp, t]. Knowledge of T[tp, t] contains all information on the time evolution of 
correlation functions. In particular, the mean value < Xm{t) > corresponds to the 
(time-dependent) value of tp m for which T has its minimum. We should emphasize 
that in contrast to the equilibrium effective action the non-equilibrium effective 
action F[ip, t] depends on the specific probability distribution for the initial values 
X° at t = to, as specified by So[x ] or, equivalently, by r[<p, to]. Nevertheless, some 
features, as the approach to a possible equilibrium state for t — > oo, can become 
independent of the initial conditions. 

The time dependence of F[<p, t] is specified by an exact evolution equation (d t is 
the time derivative at fixed ip) 

m<p\ = -P^m (6) 



fmniri2,— 713 (J^ )nin3 "3 (J^ )ni,— 

)nipA ,n2P2K >pm3 dip* pi dip* P2 dip p J 
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+ f mn l ...n k _ 1 ,-n k ( l Pn 1 + ^) mpi o * ) 

fc=4 ^Pi 

-(vn fc _ a + (r (2) ).- 1 _ 2Pfc _ 2 ^-)(^ fc _ 1 < + (r (2) )i in J (7) 

Here the matrix T^ 2 ^ [</?] is given by the second derivatives 

(r«) m „ = (8) 

and its inverse (T^) -1 denotes the effective (^-dependent) propagator. Eq. (6) is 
the central equation of this letter and we observe that it simplifies considerably if 
the equation of motion contains only up to quadratic (k = 2) or cubic (k = 3) terms. 
In the limit of infinitely many degrees of freedom Xm(t) = x( x it)i ^ * s a functional 
differential equation and d/dip m may be replaced by 5/5ip(x). Eq. (6) follows from 
the evolution equation for Z[j,t] (with d t at fixed j here) 

dtZ[j) = f m F m [^]Z[j] (9) 

and the identity d t T\ v = —d t \nZ\j. In turn, eq. (9) can be proven by taking a time 
derivative of eq. (4), inserting the equation of motion (1) under the (functional) 
integral and expressing the powers of Xn t (t) by derivatives with respect to j*. . 

The evolution equation for the mean value <p m (t) =< Xm(t) > obtains from the 
minimum condition (dF / dip m )(<f(t)) = which is valid for all t 

j t <p m (t) = F m m)) (io) 

Comparing the "macroscopic" equation of motion (10) which includes the fluctua- 
tions with the microscopic one (1) for a particular configuration, we see that the influ- 
ence of the fluctuations appears through the terms in F m (7) which involve the prop- 
agator (I^ 2 )) -1 and those involving higher (^-derivatives of T. The mean values of the 
connected two-point correlation function < Xm(t)x n (t) > — < Xm(t) >< > 
= (r^)^(^(i)) evolves according to 

-^Pmn = (</?(£)) = A pm P pn + P mp A pn 

In the limit where dF m [ip]/dip n is independent of ip (in particular for a linear equation 
of motion where A mn = —fj^l) one can solve eq. (6) explicitly. The time evolution 
of all 1PI n-point functions is given by appropriate contractions with A similar to 
eq. (11) and therefore determined by the eigenvalues of A. For small nonlinearities 
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one can perturbatively expand around this solution. Of course, the magnitude of 
the nonlinear effects is not only determined by the values of the coefficients but 
also by the size of the fluctuations. 

Let us next come to the interesting question if the evolution equation (6) has fixed 
points r*[</?] which are approached for t — > oo. Such fixed points could correspond 
to possible equilibrium distributions to which the system is attracted. The existence 
of fixed points (solutions of d t r[<£>] = 0) obviously depends on the precise form of the 
equation of motion (1). We concentrate in the following on conservative equations 
of motion 

Xu-B st —, X*s--B stMt (12) 

which can be derived from some generalized potential U[x] which is conserved if U 
does not depend explicitly on time (dll/dt = d t U\ x ). For conserved U one can show 
that 

Z*\j,0\ = J D X exp(-(3U[ X ] + f mXm ) (13) 

is a fixed point of eq. (9). (Here m = (1, s),(2, s) is a collective index.) The 
Legendre transform (5) r* is then a fixed point for the evolution equation (6) 
for arbitrary values of (3. The proof relies on the invariance of the measure D\ under 
variable shifts x' m = Xm + e m wich constant (infinitestimal) e m . Expanding 

OO 

U b(\ = E TT< ) ...n ft Xn 1 Xn 2 ...Xn fc (14) 
fe=0 K - 

the shift invariance implies the Schwinger-Dyson equation ||] 

P £ = jm (15) 

The coefficients f( k ' in the equation of motion (1) are directly related to -u( fc+1 ) and 
the right-hand side of eq. (9) vanishes due to the relative minus sign in (12). 

The symplectic structure of eq. (12) is typical for second-order equations of 
motions 

Xm = ( 16 ) 

Identifying xi with x an d X2 with it = x the equation of motion takes the form (12) 
with B st = 6 st and 

U = V[ X ] + i<7r m (17) 

We interpret U as the total energy of the system and (3 as the inverse temperature. 
The "kinetic part" of which involves the sources for 7r m factorizes as a trivial 
Gaussian. The remaining part of is the partition function for a system x-m with 
(potential) energy V[x] in thermodynamic equilibrium at temperature 1//3. We note 
that the concept of temperature arises here without any notion of interaction with 
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some external heat bath! Since U is conserved and < U >= — din Z^j = 0]/d/3 the 
temperature is uniquely fixed by the initial value < U > [t ] if Z*[0, (5] is monotonic 
in the appropriate range. 

If the thermodynamic equilibrium partition function were the only fixed point 
of the system and attractive for t — > oo, we could be sure that initial conditions 
not too far from equilibrium would "thermalize" , i.e. T[tp, t — > oo] = r* [tp]. The 
approach to the equilibrium fixed point may, however, be obstructed by the existence 
of conserved quantities A { besides U. The fixed point can be reached only if the 
initial values of all Ai coincide with the equilibrium values. The generic situation 
may be similar to quantum mechanics where the trace of arbitrary powers of the 
density matrix is conserved. In what sense (reduction of effective degrees of freedom, 
time averaging) an equilibrium state can be reached under such circumstances is an 
interesting and subtle question. 

Let us finally turn to more practical applications of eq. (6) to non-equilibrium 
situations. Exact solutions are out of reach and one has to ressort to approximations 
by truncating the most general form of F[ip, t] to an object of manageable size. The 
truncation has to be adapted to the particular system and the success of the method 
in the nonlinear non-equilibrium domain will depend on the ability to conceive a 
realistic truncation. As an example we consider the O(N) symmetric y2 4 -theory in 
d dimensions, with equation of motion (i = l..d,a = l...N,x a (x,t) real) 

X a {x, t) = didiXa{x, t) - m 2 Xa{x, t) - ^\Xb{x, t)xb{x, t)xa(x, t) (18) 

For A — > oo with fixed negative m 2 /\ this corresponds to Ising or Heisenberg models, 
and for d = 3, iV = 4 and m 2 < this is the equation used in numerical simulations of 
pion dynamics related to the investigation of a possible disordered chiral condensate. 
For the non-equilibrium effective action we use the truncation 

r>,t] = ^/^{^)^(g)^(g)+5(g)<(g)^(g) + 2C'(g)7r:(g)^(g)} 

+ J d d x{u p 2 {x) + vp(x)7r a (x)p a (x)} (19) 

with p = ^ip a <p a , ip a (q) = J d d xe mXl p> a (x) = <p*(—q). The evolution equation (6) 
describes the time dependence of the real functions A(q), B(q),C(q) which depend 
here only on the invariant and the constants u and v. Inserting (19) in (6) and 
omitting on the r.h.s. terms ~ v, one finds after some algebra 

d t A{q) = 2u; 2 C(q), d t B(q) = -2C(q) 

d t C(q) = UJ 2 q B(q) - A(q) 

d t u = 4\C(0)K 

d t v = -2u + 2\B(0)K (20) 
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with 




-G{q')G{q")G{q + q'-q")) 



+ 2 \/3J J (27r) d (27r) 
B(q)/(A(q)B(q) -C 2 (q)) 



3iV + 20 (u\ 2 f d d q d d q' 



- d G 2 {q)G{ q ')G{ q -q>) 



G(q) 



(21) 



We assume here that the system is regularized by a momentum cutoff for q 2 . For a 
given initial distribution specified by A, B, C, u and v for t = to the system (20) (or 
a simplified version with B(q) = B(0),C(q) = C(0),A(q) = A(0) + Zq 2 ) could be 
solved numerically. (Note that eq. (20) implies that A(q)B(q) — C 2 (q) is a conserved 
quantity for every q.) Results may be compared with direct numerical simulations 
which solve eq. (18) numerically and average over initial conditions. In particular, 
the expectation value of a space-independent field obeys (p = |<p a <Pa) 



or (p a = 0. It is easy to verify that the system (20) has the equilibrium fixed 
point (13) with C*(q) = 0, = 0,B*(q) = 1 and A*(q) = u} 2 ,u* = XK. The two 
last conditions are the (classical!) Schwinger- Dyson equations for the propagator 
and the quartic scalar coupling in the appropriate truncation. In thermodynamic 
equilibrium the system exhibits spontaneous symmetry breaking for A*(0) < 0. 
The contributions involving C an v on the r.h.s. of eq. (22) describe the non- 
equilibrium deviation of p from its equilibrium value p* = — A*(0)/w*. Not too 
far from equilibrium B(0) is positive and the non-equilibrium effective mass term 
A(0) — C 2 (0)/-B(0) is smaller than in equilibrium, thus enhancing the tendency to 
spontaneous symmetry breaking. 

We finally remark that the present formalism can be extended to stochastic equa- 
tions of motion and the description of correlation functions at unequal time. This 
will permit to establish the contact to path integral formulations @] and quantum 
field theory. 
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